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We present an experimental study of the directional-solidification patterns of a nematic - smectic 
B front. The chosen system is C4H9 — (CeHio^CN (in short, CCH4) in 12^im-thick samples, and 
in the planar configuration (director parallel to the plane of the sample). The nematic - smectic 
B interface presents a facet in one direction -the direction parallel to the smectic layers- and is 
otherwise rough, and devoid of forbidden directions. We measure the Mullins-Sekerka instability 
threshold and establish the morphology diagram of the system as a function of the solidification rate 
V and the angle 80 between the facet and the isotherms. We focus on the phenomena occurring 
immediately above the instability threshold when 8q is neither very small nor close to 90°. Under 
these conditions we observe drifting shallow cells and a new type of solitary wave, called "faceton", 
which consists essentially of an isolated macroscopic facet traveling laterally at such a velocity that 
its growth rate with respect to the liquid is small. Facetons may propagate either in a stationary, 
or an oscillatory way. The detailed study of their dynamics casts light on the microscopic growth 
mechanisms of the facets in this system. 
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I. INTRODUCTION 

A crystal growing from an undercooled melt rejects 
heat and chemical species, which must diffuse away in 
the liquid for the process to continue. The thus gener- 
ated thermal and solutal gradients tend to destabilize the 
advancing solid-liquid interface. This effect is counterbal- 
anced by the surface tension and the so-called intcrfacial 
kinetics, which tends to slow down the progression of the 
interface, and hence stabilize it. As a result of the com- 
petition between these conflicting factors, solidification 
fronts may assume a large variety of nonlinear patterns, 
the characteristics of which depend on the control pa- 
rameters, and the initial and boundary conditions of the 
process. 

The study of solidification patterns has been an active 
field of research for several decades jfl-[|. Most of the 
existing studies are devoted to fully non-faceted systems. 
In such systems, the surface tension 7 and the kinetic co- 
efficient (3 (defined as the ratio of the kinetic undercooling 
to the growth velocity) are non-singular functions of the 
orientation of the interface with respect to the crystal lat- 
tice. On a molecular scale, this corresponds to the fact 
that the interface is rough in all orientations. Familiar 
aspects of the dynamics of fully non-faceted systems in 
directional solidification (i.e., when the system is pulled 
at a constant velocity V toward the cold side of an ap- 
plied unidirectional thermal gradient G; see Fig. [l]) are 
the existence of a stable planar front at low values of 
V, the primary cellular (or Mullins-Sekerka) instability 
occurring at a threshold velocity V c , the quasi periodic 
arrays of rounded cells at V slightly above V c , and of 



dendrites at V much higher than V c . Many dynamical 
features of these patterns (e.g., stability limits, modes 
of instability) are not yet fully understood, but some of 
their fundamental properties are now clear, among which 
the crucial role played by interfacial anisotropy fj])(i]|| ■ In 
fact, a certain minimum degree of interfacial anisotropy 
is a necessary condition for cellular and dendritic arrays 
to be stable, or even exist. In thin samples -i.e., quasi 
bidimensional (2D) systems- , 7 and (3 are functions of a 
single variable, say, the tilt angle 9 of the interface with 
respect to the isotherms. The functions 7(6*) and (3(9), 
and thus the solidification patterns, depend on the ori- 
entation of the crystal with respect to the solidification 
setup 
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FIG. 1. Sketch 
of a thin-sample directional-solidification setup, z: axis 
of the thermal gradient, x: axis parallel to the isotherms. 
V: pulling velocity. After a transient, the front advances 
(in average) at the imposed velocity V with respect to the 
liquid, and thus remains essentially immobile in the labora- 
tory reference frame. It can then be continuously observed 
with an optical microscope. 

In contrast with the case of fully non-faceted systems, 
little is yet known about the directional-solidification dy- 
namics of faceted crystals. The few existing experimen- 
tal studies on this subject first of all show that a dis- 
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tinction must be made between fully and partly faceted 
systems |6||^-|l2]]. Growth facets (which most generally, 
although not necessarily, coincide with equilibrium facets 
) correspond to planes of the crystal containing sev- 
eral directions of strong binding. Fully faceted crystals 
have numerous facet directions, and their directional- 
solidification fronts consist of a succession of facets lim- 
ited by sharp edges. The dynamics of such fronts does 
not give rise to any stationary state, in general, and 
bears no obvious relation with that of non- faceted fronts. 
Partly faceted systems only have a few facet directions 
connected to one another by large rounded regions. In 
lamellar crystals, the solid-liquid interface may be rough 
in all but one direction, namely that of the molecular 
layers. In this case, when the tilt angle 9q of the layers 
with respect to the isotherms is large, the dynamics of 
the front must obviously be that of a non-faceted crystal 
as long as the deformation of the front remains small, 
that is, below V c and in a small range of V above V c . 
Facets only appear at higher V when the deformation 
of the interface is large. A relatively smooth transition 
from the non-faceted to a partly faceted dynamics may 
then be observed. This is the experimental configuration 
considered in this study. 

In this paper, we study the directional-solidification 
dynamics of the front associated to the nematic-smectic 
B transition of the liquid crystal C4H9 — (CeHio^CN (in 
short, CCH4). A long-range order exists in the direc- 
tion perpendicular to the molecular layers in the smectic 
B phase, so that this phase actually is a lamellar crys- 
tal. Previous free-growth studies have indeed shown that 
the nematic-smectic B fronts of the n — 3, 4, 5 members 
of the series CCHn (where n stands for the number of 
carbon atoms in the aliphatic chain) have a single facet 
direction parallel to the molecular layers of the smec- 
tic phase, and are rough in all other directions [O-M. 
Moreover, they have no unstable orientations in a di- 
rection perpendicular to the molecular layers, contrary 
to the Smectic A-Smectic B fronts previously studied in 
directional solidification by Oswald et al [|6| JTl|Jl2]] . The 
present study is performed in thin (12 /xm-thick) samples 
and in the planar configuration (director parallel to the 
plane of the sample), in order for the front -including 
the facets, if any- to remain perpendicular to the sample 
plane. Practically, the system is thus a 2D one. 

We shall mostly focus on a type of solitary wave ap- 
pearing near the Mullins-Sekerka threshold, called "face- 
ton" because it contains a single small facet traveling 
along the front at such a velocity that the normal growth 
rate of the facet, i.e., its growth rate with respect to the 
liquid, is generally much smaller than V. Such a phe- 
nomenon, which has never been observed before, to our 
best knowledge, is obviously highly specific to faceted di- 
rectional solidification, and therefore particularly inter- 
esting from our present viewpoint. A preliminary com- 
ment about the nematic-smectic B facets in the CCHn 



series is in order. The growth rate of a facet is con- 
trolled by the dynamics of the molecular steps flowing 
along it. Therefore it crucially depends on whether, or 
not, the facet contains, or is connected with step sources 
jD|[L7|. When no step source is available, the facet grows 
through nucleation and spreading of new terraces (sur- 
face nucleation), which is a very slow process at low un- 
dercooling. In fact, the growth rate of a perfect facet 
is totally negligible when the undercooling is lower than 
some finite value. Such a behavior ("blocked" facets at 
low undercoolings) has clearly been observed during the 
solidification of many, but not all the studied faceted sys- 
tems. What concerns us here is that it was not observed 
during the free growth of the smectic B phase of CCH3, 
despite the strongly faceted aspect of the growing crys- 
tals 0-|l6| . Numerical simulations in which a cusp-like 
minimum of "f{9) but no anisotropy of (3 was taken into 
account satisfactorily reproduced the observed growth 
shapes. Thus the observation of a facet on a macroscopic 
scale would not necessarily mean the presence of a sin- 
gularity in (3. In order to clarify this point in the case 
of CCH4, we report, in Section III, preliminary observa- 
tions in free growth showing that the nematic-smectic B 
facet of CCH4 is capable of remaining immobile at un- 
dercoolings lower than O.IK. Thus, in CCH4 at least, 
the nematic-smectic B fronts can form growth facets. 



II. EXPERIMENTAL 

The relevant material parameters of the liquid crys- 
tal CCH4 (MERCK IS-0558) can be found in ref. |n|. 
The residual impurities, the chemical nature of which is 
unknown, were characterized as regards solidification by 
the usual methods (see below). We found that the im- 
purity content at the outset of the experiments was re- 
producible, but slowly increased during the experiments, 
indicating that the product was undergoing a decompo- 
sition in the nematic phase, as previously noticed and 
analyzed for the case of CCH3 Jll| . The nematic - smec- 
tic B transition temperature TVs wa s generally of about 
53.1K in fresh samples. Fig. || shows TV 5 measured as 
a function of time in one sample. It can be seen that 
the decomposition rate is sufficiently slow not to severely 
perturb a solidification run, but sufficiently rapid to pre- 
vent us to carry out several successive runs with the same 
sample. Outgasing the as received product resulted in a 
significant slowing down of the decomposition process. 

We have studied the crystal structure of the smectic 
B phase of CCH4 by low-angle X ray diffraction Jl9| . 
As expected, this phase is basically an AB-type stacking 
of hexagonal layers. The parameters are approximately 
a = 5.9A and c = 29A, which is in accordance with the 
data available for the other members of the homologue 
series pfl]. The hexagonal layers however appear to be 



2 



slightly distorted, which may entail the existence of su- 
perstructures in the layers. 

A schematic view of a thin-sample directional- 
solidification experiment is shown in Fig. [P. A detailed 
description of our setup is given elsewhere . In this 
study, the samples were made of two parallel glass plates 
separated by 12 /xm-thick plastic spacers. Their useful 
width was of 9 mm and their length of 60 mm. They 
were filled under an argon atmosphere at a temperature 
higher than T/vSj an( i then cooled down to room tem- 
perature. Numerous smectic crystals appeared by het- 
erogeneous nucleation during cooling. The samples were 
placed in the thermal gradient, and a smectic B crystal 
of known orientation (determined through the observed 
value of 0o ) was selected by a method to be explained 
shortly. The sample was annealed at rest for about 30 
minutes in order to homogenize the concentration in the 
liquid. V was then switched to a chosen value, left for 
given time at this value, and then increased step by step. 
The temperature gradient at the growth front was of 53 
Kcm ■ , unless otherwise mentioned. The pulling veloc- 
ity was in the range of 0.3 - 30 fims . The observations 
were made with a polarizing microscope (Leica) equipped 
with a CCD camera. The video signal was analyzed with 
digital image processing. 
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FIG. 2. The nematic - smectic B equilibrium temper- 
ature in a CCH4 sample as a function of time. Tns was 
measured by controlling the temperature of a free-growth 
stage in order to keep a small smectic B crystal in quasi 
equilibrium with the nematic. The relatively low initial 
value of Tns indicates that the sample was rather impure 
at the outset. 

It is obviously crucial for our experiments that large 
smectic B crystals of arbitrary orientation might be se- 
lected. To this aim, we have studied the influence 
of various treatments of the inner sides of the glass 
plates. Three types of plates were used : untreated 
plates, plates covered with a mono-oriented thin film 
of poly(tetrafluorethylene) (PTFE) prepared by friction 
transfer at T « 200 °C §, or with a w lOOA-thick 
layer of Al or In deposited by oblique evaporation. In 
the nematic phase, the orientation of the director was es- 
sentially planar in all samples. The director was more or 



less aligned along the direction of friction, or deposition, 
in treated samples, but domains corresponding to small 
(a few degrees) variations in the orientation of the direc- 
tor, still existed (see Fig. ^ below). This inhomogeneity 
of the nematic phase caused but minor perturbations in 
our experiments, since the phenomena of interest turned 
out to be essentially independent of the orientation of 
the nematic director. In all samples, the smectic B phase 
had a planar orientation, but was divided into different 
crystals (or grains) corresponding to a different value of 
0o . The surface treatment gave a pronounced preferential 
distribution of 0o among the various grains, facilitating 
the selection of the desired value of 0o- The size of the 
selected smectic B grain was increased by a method con- 
sisting of forcing the crystal to grow through a funnel 
shaped obstacle fjj. By this method, smectic grains of 
a millimetric width, and arbitrary values of 0o were ob- 
tained. 



III. CHARACTERIZATION OF THE SYSTEM 
A. Free growth at small undercoolings 

The observations reported in this Section were per- 
formed with a free-growth setup similar to the one de- 
scribed in ref . |2^] , in which the changes in the undercool- 
ing are produced, via the Clausius-Clapeyron effect, by 
a sudden pressure change at constant temperature, and 
are therefore quasi instantaneous p2] , p3| . The samples 
were the same as those used in directional solidification. 
At the beginning of the experiments, the samples were 
heated step by step until only one small smectic B crys- 
tal was left in the nematic. The sample was maintained 
at constant temperature until the changes in the shape 
of the crystal became very slow (this took about 20 min- 
utes). Admittedly, this shape is not the exact equilibrium 
shape of the crystal, but it exhibits clear reproducible 
features, namely, long facets parallel to the smectic lay- 
ers and rounded ends in the perpendicular direction ( 
Fig. ^| a.), which is enough for our present purpose. 
It should be noted that the observed near-equilibrium 
shape clearly shows the absence of a forbidden orienta- 
tion range around 89 = 90°, where 89 is the deviation of 
the interface from the direction of the molecular layers, 
but suggests that the facet might actually be limited by a 
sharp edge, i.e., the interface might be unstable at small 
values of 89. The fact that we have not observed the Her- 
ring instability [ p4|JTT| in directional solidification at the 
lowest explored value of 0o indicates that this forbidden 
orientation range is very small (< 2°), if it exists at all. 

A sudden increase of the undercooling 8T was applied 
at time t = 0, and the subsequent growth of the crystal 
recorded (Figure ^). The growth process, which is gov- 
erned by the anisotropic interfacial properties and dif- 
fusivities, is very complicated. Its study is beyond the 
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scope of this paper. Here we limit ourselves to the fol- 
lowing observation: the facets of the smectic B crystals 
remained blocked within experimental uncertainty (their 
growth rate was lower than about 0.01/ims" 1 ) at under- 
coolings lower than 0.1 K (Fig. |^). At higher under- 
coolings, they generally grow at a measurable rate. The 
apparent threshold undercooling AT nuc i for growth by 
surface nucleation of our system is thus larger than 0.1 K 
and probably not much larger than this value. This esti- 
mate of AT nuc [ is small compared to what it is in ordi- 
nary solid-liquid systems, but this may be explained by 
the small value of 7 in our system |2^,|26| . It is also possi- 
ble that, in our thin samples, surface nucleation is in fact 
heterogeneous, i.e., takes place preferentially along the 
line of contact with the glass plates. The nucleation rate 
would then depend on the treatment of the glass plates. 




FIG. 3. Free growth. Successive snapshots of a smec- 
tic B crystal of CCH4 growing from the nematic phase at 
AT = 0.07^. 



B. Directional solidification: instability threshold 

The Mullins-Sekerka instability threshold was found to 
lie between approximately 2 and 3 /ims -1 in all the stud- 
ied (fresh) samples. No influence of the orientation of the 
smectic, or the nematic was observed within experimen- 
tal uncertainty. However, it should be noted that this 
uncertainty was large (« 1 ^ims~ l ) for the reason to be 
explained presently. 

Figure [| shows a sample at rest, and pulled at a rate 
lower than V c . Two isothermal fronts are visible, namely, 
a front separating the nematic (N) phase from a smectic 
B domain (Sml), and at a lower temperature, a front 
separating Sml from a second smectic B domain (Sm2). 
The nature of the transition from Sml to Sm2 is not yet 
clear. This transition was observed in most, but not all 
experiments. Observations (not reported here) incline us 
to think that Sm2 is the same phase as Sml, but with a 
different orientation, thus a different interaction energy 
with the glass plates. In any case, we need not take into 
account the Sml-Sm2 front here since this front, when 
present, does not perturb the dynamics of the N-Sml 
front. 

It can be seen in Fig. |^ that the nematic- smectic B 
front remains planar during solidification at V < V c , ex- 



cept for small, long-wavelength distortions due to the 
presence of domains in the nematic phase. These dis- 
tortions are larger during solidification than at rest, and 
undergo sudden changes each time the front leaves a ne- 
matic domain for another. This phenomenon has thus 
an equilibrium as well as a kinetic origin. In our ex- 
periments, it plays the role of a relatively strong long- 
wavelength, low-frequency noise, which blurs some of the 
morphological-transition thresholds of the system. This 
is the main origin of the aforementioned large uncertainty 
on the measured values of V c . However, we can state with 
certainty that V c was higher than 2 /ims -1 since the dis- 
tortions caused by nematic domains, or any other source 
of perturbation (e.g., dust particles) did not amplify be- 
low this velocity. 
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FIG. 4. Directional solidification (in this, and all the fol- 
lowing micrographs, growth is upwards). N: nematic. Sml: 
smectic B. Sm2: smectic B oriented differently from Sml. 
(a) Sample at rest (V = 0); (b) Sample in the process of 
solidification at V = 0.9/xms -1 . Note the domains in the 
nematic. Sml is a single crystal, but Sm2 a polycrystal, as 
shown by the presence of cusps on the Sml-Sm2 front. 



C. Solute redistribution transient 

When V is smaller than V c , the front reaches a station- 
ary planar state through the so-called solute redistribu- 
tion transient. A recoil curve -i.e., the curve representing 
the variation of the position (or temperature) of the pla- 
nar front as a function of time during the initial transient 
of a particular run- is reproduced in Fig. |^. It is well 
known that information about the relevant properties of 
the solute (diffusion coefficient D in the liquid, partition 
coefficient K, thermal gap AT Q ) can be gained from the 
characteristics of the transient, and the value of V c . We 
have utilized this method in order to characterize the un- 
known impurity playing the role of solute in our system. 

The threshold velocity, and the amplitude of the so- 
lute redistribution transient are given by V c « (1 + 
KD s /D)DG/AT a {D s : diffusion coefficient in the solid) 
and AT , respectively | p7fl . By fitting the recoil data us- 
ing the Warren-Langer approximate theory [E8| (Fig. ||) 
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and assuming V c = 2.5^,ms~ 1 and KD S /D << 1, we 
obtained K = 0.12. This gives D = 80 fxm 2 s~ 1 and 
AT = 0.2 iC. These data give us no information about 
D s , but there is good reason to believe that our system 
is a two-sided one -i.e., that D s is not much smaller than 
D g§. 

o I . 1 1 1 . 1 o 



It can be seen that the sequence of morphologies ob- 
served as a function of V for a fixed value of 6q is the 
same for all values of 9 , except for those close to 0° 
(facets parallel to the growth front) or 90° (facets per- 
pendicular to the growth front). This generic sequence is 
illustrated in Fig. ^. 
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onset of the recoil is an instrumental effect. 



FIG. 5. Recoil curve at V 
as in Fig. ^| Continuous line: 



IV. RESULTS 



A. Morphology diagram 



A diagram displaying the observed morphologies as a 
function of the pulling velocity and the orientation of the 
smectic B crystal is shown in Fig. ^. 
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FIG. 6. Morphology diagram. Measurement points: 
waves and facetons (o), facetons and unstationary faceted 
fingers (□), unstationary faceted fingers (A), stationary 
faceted fingers (•) and unstable facets (x). Heavy dashed 
line: Mullins-Sekerka instability threshold. Inset micro- 
graphs : see Fig. 0. 




FIG. 7. The different growth morphologies observed 
as a function of V for 8q — 25°. (a) Planar front; 
(b) Drifting shallow cells; (c) Drifting faceton (stationary 
mode); (d) Drifting facetons (oscillatory mode) at differ- 
ent stages of their oscillation cycle; see Fig. [l8| below; (e) 
Non-stationary array of faceted fingers ; (f) Stationary ar- 
ray of faceted fingers. 



Small-amplitude, nearly sinusoidal traveling waves ap- 
pear near the instability threshold (Fig. 0.b), in accor- 
dance with previous observations in two-sided anisotropic 
systems ||. Such weakly nonlinear waves are commonly 
called " shallow cells" . 

We observed drifting shallow cells in a broad range of 
V around the threshold (1 /jms -1 < V < 8 /ims -1 ). In 
the same range of V, we also observed "facetons" (Fig. 

c). These solitary waves can propagate in a stationary 
or an oscillatory way. They appear when the amplitude 
of the cells is so large that the tilt angle of the front locally 
reaches the value Qq corresponding to the facets. Most 
generally, this occurs under the effect of perturbations 
due to the nematic domains. The frequency of creation 
of facetons, and thus their average number by unit length 
of the front increases as V increases. When the average 
spacing of the facetons becomes smaller than their width 
(ss 200/xm), they cease to behave as non- interacting ob- 
jects. In fact, they disappear altogether, giving way to 
arrays of much narrower objects, called faceted fingers 
(Fig. |7[e). This occurs at about 8/zms _1 . However, this 
transition is strongly noise dependent, and thus relatively 
ill-defined from an experimental viewpoint. Shallow cells 
and facetons are studied in detail in the next section. 
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FIG. 8. Transition from isolated facetons to faceted fin- 
gers for So = —70° in an aged sample, (a) V = 3.1 fims' 1 
(snapshot of the front); (b) Corresponding spatiotem- 
poral diagram (time series of the intensity distribu- 
tion along a line located 20 /im below the front); (c) 
V = 6.5 fj,ms~ ; (d) Corresponding spatiotemporal di- 
agram; (e) V = 13.5 (jros -1 . Note that another grain 
(#o = 73°) appears in the leftmost part of the figure; (f) 
Corresponding spatiotemporal diagram. 

The arrays of faceted fingers, which are observed above 
8 /ims" 1 exhibit a relatively sharp transition from an un- 
stationary (Fig. [?].e and Fig. |J) to a stationary dynamics 
as V increases (Figs. [j].f and the dispersion appear- 
ing in Fig. ^ is mostly due to the aging of the samples) . 
The spatiotemporal diagrams of the unstationary arrays 
shown in Fig. |^ reveal the transitory or local existence 
of well-defined oscillatory modes. These modes become 
more and more apparent as V increases because the os- 
cillation period T osc is a rapidly decreasing function of V 



(Fig. |T^.a). This strongly suggests the existence of a ho- 
mogeneous oscillatory bifurcation of the high-F station- 
ary patterns as V decreases within some narrow range of 
spacing. 
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FIG. 9. Stationary array of faceted fingers at 
V = 13.5 /^ms _1 and 6o = 24° (a.) Snapshot of the front; 
(b.) Spatiotemporal diagram (piling- up of skeletonized im- 
ages of the growth front). 
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FIG. 10. Oscillation period (a) as a function of V for 
8o = 56° (b) as a function of |#o| for three values of V. 
The leftmost point in a) corresponds to an isolated oscilla- 
tory faceton. 

We now turn to the particular orientations correspond- 
ing to the bounds of the scanned interval of 6$. When 
6q = 90°, the system is reflection symmetric. Shallow 
cells no longer drift, and facetons cease to exist. The 
shallow cells break up into narrow faceted fingers as V 
is increased above threshold (Fig. [Ilia). The widest 
faceted fingers, which are the majority ones, are not re- 
flection symmetric, whereas the narrowest ones are reflec- 
tion symmetric. The two opposite but equivalent direc- 
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tions of symmetry breaking are equally populated. The 
resulting arrays were non stationary even at the highest 
explored values of V (Fig. |Xl|.b). This is very differ- 
ent from what was observed by Oswald et al in smectic 
A-smectic B fronts for a similar orientation of the facet 
|@|. In that system, because of the existence of forbidden 
directions, the finger tips exhibited pointed triangular 
shapes, and formed stationary arrays. 



At higher V, they evolved as illustrated in Figs. [TJ] and 
111 
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FIG. 11. (a) Array of symmetry-broken faceted fingers 
at #o = 90° and V = 10 urns' 1 ; (b.) Corresponding spa- 
tiotemporal diagram; (c) Instability of a facet at 9o = 2° 
and V = 10 /ims -1 . 

When #o is sufficiently close to zero, the growth front 
of smectic B grains is entirely occupied by a facet at any 
value of V. This can be considered as a finite-size ef- 
fect resulting from the following fact: facets are always 
present in the grooves attached to grain boundaries for 
whatever values of 9o and V; the stationary size of these 
facets is more or less proportional to l/(tan\9o\); they 
thus occupy the whole grain when \9q\ is lower than a 
certain value, which is of about 2° for a grain size of 
500fim. At sufficiently high V, these long facets break 
up through the mechanism illustrated in Fig. |ll].c. It 
is not necessary to repeat here the description of this 
process, which has been presented by other authors [ jl2| . 
We simply note that, in our fresh samples, this instabil- 
ity was observed to result from the occasional collisions 
of the front with defects (domain walls, dust particles) 
present in the nematic. In the less pure samples, it was 
superseded by another well-known process, namely, the 
nucleation of crystals in the undercooled melt ahead of 
the front Both mechanisms give rise to more or less 
permanently cyclic growth regimes. 



B. Near-threshold patterns 

1. Drifting shallow cells 

Most generally shallow cells appeared in the form of a 
noise-induced wave packet. A spontaneous homogeneous 
growth of the cells was never observed with certainty. We 
attribute this fact to the interplay between shallow cells 
and facetons (see below). At, or below 2fims~ 1 , noise- 
induced wave packets systematically disappeared when 
the source of noise disappeared, as already mentioned. 



120 s 



100 urn 



60s 




Spatiotemporal diagram of a drifting wave 
packet; V = 3.1 /ims' 1 , 8o = 25°. 

A careful analysis of the spatiotemporal diagram of 
Fig. |l2] has shown that (i) the cells are initially sinu- 
soidal; (ii) they grow in amplitude with a uniform am- 
plification rate of « 0.002s -1 ; (iii) the amplest cells are 
no longer sinusoidal at the end of the time sequence, (iv) 
the spacing A and the drift velocity Vd are uniform in 
space and constant in time. Vd is thus amplitude in- 
dependent. This is in keeping with the idea that this 
sequence is the initial stage of the usual amplification 
process leading from a linearly unstable state to a sta- 
tionary weakly nonlinear regime. The final regime was 
not observed because the process was interrupted by an 
external perturbation giving rise to a faceton. 

The traces on the lefthand side of Figure [l3| are the tra- 
jectories of three oscillatory facetons. These objects are 
studied below. For now, the point of interest is that the 
rearmost faceton leaves behind a region of the front which 
is free of detectable shallow cells (see also Figs. |l^ and [l8| 
below) . The cells reappear at w 200 fim from the faceton, 
and then amplify following a process entirely similar to 
the above one, except for two points: (i) in the present 
case, the amplification rate (« 0.02s -1 ) is much larger 
than in the preceding case, since V is higher, and (ii) a 
stationary regime of nonlinear shallow cells is reached. 
This confirms clearly, although only scmi-quantitatively, 
that the system admits stationary weakly nonlinear cel- 
lular states within a measurable range of V above V c . 
These states are metastable with respect to the forma- 
tion of facetons. Also we note that the direction of drift 
of the cells is opposite to that of facetons. This is some- 
what of a surprise since, in other systems, shallow cells 
and facets have been found to drift in the same direction 
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100 ran 




FIG. 13. Spatiotemporal diagram. The three traces on 
the lefthand side are the trajectories of oscillatory face- 
tons drifting leftwards. Note the disappearance of the cells 
(which drift rightwards) in the wake of the rearmost face- 
ton. A temporary exception to this rule is visible near the 
end of the recording, when the faceton emits a packet of 3 
or 4 cells. This exception is only apparent, however, since 
this occurs during a period of time when the faceton no 
longer exists (it is drifting rightwards). V = 6.5 ferns' 1 , 
6o = 55°, recording time = 250 s. 

The measured values of \Vi\ and A are plotted in Fig. 
|l4] as a function of 0q for a given value of V. The data 
are compatible with the fact that Vd(0o) must go to zero 
at 6{) = 0° and 90° for symmetry reasons. The maximum 
is at about 70°, and corresponds to a relatively large 
value of Vd/V, indicating that the system is strongly 
anisotropic even in the orientation range in which the 



interface is rough. 
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FIG. 14. Drift velocity (a) and wavelength (b) of the 
cells as a function of the tilt angle of the facet at 
V — 6.5/ims _1 . 

We have noted above that Vd seems to be independent 
of the amplitude of the cells. It is thus legitimate to ad- 
mit (but not certain) that the measured value of Vd is the 
same as in the linear regime. We have performed a linear 
stability analysis of the planar front of a two-sided sys- 
tem taking into account the anisotropy of the diffusion in 
the two bulk phases (nematic and smectic B), and that of 
the linear kinetic coefficient (3 (the anisotropy of 7 does 
not come into play in a linear calculation Q). We have 
solved the dispersion equation numerically under various 
assumptions concerning the orientation dependences of 
D, D s , and /3, which are not known. Qualitatively, the 
results can be summed up as follows ^7|. We find that 
the observed sign and absolute value of Vd(0o) could be 
ascribed to diffusion anisotropy only if, in the smectic B 
phase, the impurities diffused much faster through the 
smectic layers than parallel to them, which is very un- 
likely to be true. Thus the observed drift of the shallow 
cells is most probably due to kinetic anisotropy. In such a 
case, the sign of Vd is given by —d(3/d0 In conclusion, 
the observed direction of drift of the shallow cells (if it is 
really the same as in the linear regime) indicates that, in 
our system, (3 increases as 80 increases. This result poses 
no particular problem except for the vicinal domain, in 
which (3 is expected to be more or less proportional to 
the reciprocal of the step density, and hence to the re- 
ciprocal of 80 p9] . The cross-over from the vicinal to 
the rough domains as 80 increases should thus manifest 
itself through a change in the sign of Vd- It is tempting 
to assume that this cross-over corresponds to the zero 
of Vd{Oo) which perhaps appears near 12° in Fig. p"4|.a. 
However, the observation of macroscopic facets drifting 
in the same direction as the shallow cells disproves this 
assumption, and indicates that the vicinal domain is ac- 
tually very narrow in our system (see below). 
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2. Stationary facetons 

The spatiotemporal diagram of a stationary faceton is 
shown in Fig. Ea. Clearly, a faceton is a solitary wave 
consisting of a macroscopic facet and a broad rounded 
finger separated from each other by a very thin liquid 
groove. The regularity of the spatiotemporal diagram 
shows that facetons, once formed, are quite stable. In 
particular, they absorb the shallow cells which they may 
encounter ahead of themselves without being modified, 
and seem to be insensitive to the perturbations caused 
by the nematic domains. The depth of the facet -i.e., the 
distance Azf between the two edges of the facet along the 
z axis- is difficult to measure with accuracy because the 
lower edge, located near the bottom of the groove, is gen- 
erally not resolved. However, it is certain that Azf is in 
the 30 — 50fj,m range (the difference of temperature ATf 
between the two edges is thus in the range 0.15 — 0.25 K), 
and decreases as 0$ increases. The upper edge of the facet 
corresponds to a small pointed maximum of the front 
shape, but it is not possible to decide whether, or not, 
this edge is sharp on a molecular scale. The width of the 
rounded finger -i.e., the extension of the deformed region 
of the front behind the finger tip- is of about 200 jim. 
As mentioned, shallow cells do not develop in this region 
of the front. The trajectory of the faceton makes a small 
angle with the direction of the macroscopic facet, indicat- 
ing that the normal growth rate of the facet is small but 
finite. Thus the facet is not blocked, and the question 
arises as to its microscopic growth mechanisms. 



the dominant mechanism. Moreover, the fact that both 
V n and Azf are decreasing functions of 9q suggests that 
V n is essentially determined by events occurring near the 
lower edge of the facet. One may imagine either that 
surface nucleation takes place at a relatively high rate at 
this point, or that the facet is supplied with steps coming 
from the bottom of the groove where the interface is nec- 
essarily rough. In both cases, V n would be very sensitive 
to the details of the conformation of the interface in this 
region. These details may depend on the treatment of the 
glass plates, which could explain the dispersion between 
values measured in different samples. 



OV=3.1|lm/s 
□ V=6.5|lm/s 
OV=10nm/s 
_V=13.5|im/s 



1 * 
d 

> 2 



A 

□ 



□ QD a 



1 m 



20 40 60 80 

e (deg) 

FIG. 16. Normal growth velocity of facets belonging to 
facetons or arrays of faceted fingers as a function of the 
tilt angle of the facet for the indicated values of the pulling 
velocity. In the case of oscillatory facetons the minimum 
value of V n has been plotted. The data point at 8o — 2° 
corresponds to the facet shown in Fig. [llj.c prior to its 
destabilization. 




FIG. 15. Stationary faceton. 8o = 25°, V = 6.5 fims -1 . 
(a.) Snapshot of the front. The faint dark line appearing 
in the solid in the continuation of the facet is a thin liquid 
groove; see Fig. [H]. (b.) Spatiotemporal diagram. The 
normal growth rate of the facet is V„ fa 0.9 fj.ms -1 



Fig. [T(j displays a large number of values of the normal 
velocity of facets V n measured in isolated facetons as well 
as in arrays of faceted fingers for various values of V and 
6q. In spite of a large dispersion of the data, it is clear 
that V n is essentially a non-zero quantity which decreases 
as 8q increases, and increases as V increases. The regu- 
larity of the stationary facetons or arrays (see Figs. [l5| 
and ^) , and the fact that V„ is very close to zero when 9q 
is large allow us to exclude screw dislocation growth as 



3. Oscillating facetons 

Fig. |l^ shows a process of formation of facetons in 
response to a perturbation. Macroscopic facets progres- 
sively develop on one side of the shallow cells as the am- 
plitude of the latter increases. These facets first drift 
with the same velocity as the shallow cells, and then 
change their direction of drift. This change is not ac- 
companied by any modification in the orientation of the 
facets within experimental uncertainty (« 0.5°). Thus 
the same macroscopic facet can be in two different mi- 
croscopic states, or growth regimes. One of these (the 
"slow" regime) is that of the stationary state, discussed 
in the preceding section, while the other (the "rapid" 
regime) corresponds to a rough interface. As announced, 
we are thus led to assume that the cross-over from vici- 
nal to rough interfaces occurs at values of 59 lower than 

0.5° in our system. This is indeed surprising since this 
disorientation corresponds to a very low density of steps 
(less than 1 per fim), but not impossible. We also note 
that the persistence of a macroscopic facet while the in- 
terface is rough on a microscopic scale is explainable by 
the sole singularity of the 7 -plot Jig , 
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FIG. 17. Facetons appearing in response to a perturba- 
tion. Spatiotemporal diagram. V = 6.5 (ims -1 , do = 25°, 
recording time: 60 s. Note the opposite signs of the drift 
velocities of the cells and the facets. 

The existence of two different growth regimes of a 
macroscopic facet is confirmed by the fact that face- 
tons most often adopt an oscillatory mode of propagation 
(Fig. [l8]). Obviously, this oscillation consists of a more 
or less ample cycle between the aforementioned rapid and 
slow regimes. The conditions under which facetons are 
stationary, or oscillatory could not be determined. In 
fact, stationary facetons were observed much less fre- 
quently than, and always in coexistence with oscillating 
facetons. Moreover, some oscillating facetons were regu- 
lar (Fig. |l^), but most of them were irregular (Figs. [l3| 
or |l^) . It is possible that the system intrinsically admits 
stationary, periodic, and more or less chaotic facetons. 
However, the following explanation is also possible. 




FIG. 18. Oscillating faceton. 9o = 42°, V = 6.5 /ims -1 . 
(a.) Snapshots of the front at different stages of an oscilla- 
tion period, (b.) Spatiotemporal diagram. 

A careful inspection of Fig. [l8| reveals that the tran- 
sition of the oscillating facetons from a slow to a rapid 
regime corresponds to a sudden pinching-off of the liq- 
uid groove, whereas the reverse transition from a rapid 
to a slow regime consists of a progressive deepening of 
the groove. If we focus on the sole groove, this behav- 
ior is strongly reminiscent of the periodic pinching-off 
(called cusp instability) of the intercell grooves in non- 
faceted cellular fronts Q. This instability, we recall, is 
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most probably of a capillary origin (Rayleigh instability) 
PH , and very sensitive to the lattice defects which, in the 
non-faceted systems, are often attached to the groove -in 
fact, the grooves to which subboundaries (low-angle grain 
boundaries) are attached are not subject to the cusp in- 
stability ]32] . If, by analogy, we assume that the intercell 
groove of facetons, like that of nonfaceted cells, is intrin- 
sically subject to an oscillatory Rayleigh instability, we 
are led to the conclusion that the transition of the facet 
from a slow to a rapid regime is a secondary effect due 
to changes occurring in the configuration of the interface 
near the lower edge of the facet. The presence of lattice 
defects (e.g., subboundaries) emerging into the liquid at 
the bottom of the groove may hinder these changes, sup- 
pressing the oscillation. This would explain that facetons 
are much more often oscillatory than stationary. 

4- Lattice defects 

Some lattice defects (mostly, grain boundaries) can be 
detected with the optical microscope thanks to the fact 
that they create macroscopic depressions (grooves) of the 
growth front around the point at which they emerge into 
the liquid. In our system, these grooves must be partly 
faceted during solidification. We have lowered the ap- 
plied thermal gradient in some experiments in order to 
facilitate the observation of such grooves. This allowed 
us to reveal that the growth front of our system is often 
swept by very small facets, called micro-facets, certainly 
attached to lattice defects emerging into the liquid. 




FIG. 19. Spatiotemporal diagram showing shallow cells, 
oscillating faceted solitary waves, and micro-facets (arrow) , 
6 = 36°, V = 3.1 fims' 1 and G = 25 Kern' 1 . 



We observed several types of micro-facets, correspond- 
ing probably to different types of lattice defects. The 
micro-facets of the type shown in Figure |l^ were rela- 
tively easy to identify because they travel at a perfectly 
constant velocity, catching up, and running through all 



the other structures of the front, in particular, facetons. 
Their drift velocity has thus most probably the maxi- 
mum possible value, i.e., the value corresponding to to- 
tally blocked facets. They must be attached to lattice 
defects -stacking faults, or twist subboundaries- strongly 
locked onto the lamella plane of the smectic. However, 
these micro-facets seem to have but little effect on the 
dynamics of the front. They indeed provoke an instan- 
taneous slowing-down of the macroscopic "rapid" facets 
when they collide with them (see Fig. |l^), but do not 
trigger a durable transition to the slow regime. So this 
observation, whatever its intrinsic interest may be, does 
not cast light on the question of the possible role played 
by lattice defects in the dynamics of the facetons. 



V. DISCUSSION 



We have shown that the directional solidification of a 
nematic-smectic B front in the planar configuration gives 
rise to a wealth of interesting nonlinear phenomena, the 
most striking of which are the stationary or oscillatory 
"facetons" encountered in the vicinity of the Mullins- 
Sekerka threshold. These observations raise numerous 
unsolved problems concerning the microscopic growth 
mechanisms of the facets, as well as the nonlinear dynam- 
ics of the observed macroscopic patterns. An important 
question is whether these phenomena are specific of the 
nematic-smectic B fronts, or are of frequent occurrence 
in partly faceted fronts. In order to clarify this point, we 
are currently searching for similar phenomena in more 
conventional, partly faceted solidification fronts. Also 
numerical simulations based on a phase-field method are 
in progress in order to test the consistency of the nu- 
merous conjectures which we have been led to make in 
order to explain the peculiar dynamical features of the 
facetons. 
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